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Abstract 

Using the 1/Nc expansion scheme and truncating the hadronic spectrum to the lowest- 
lying resonances, we match a meromorphic approximation to the {SPP) Green function 
onto QCD by imposing the correct large-momentum falloff, both off-shell and on the 
relevant hadron mass shells. In this way we determine a number of chiral low-energy 
constants of 0{p^), in particular the ones governing SU{3) breaking in the Ki^ vector 
form factor at zero momentum transfer. The main result of our matching procedure is 
that the known loop contributions largely dominate the corrections of 0{p^) to /+(0). 
We discuss the implications of our final value f^°^ (0) = 0.984 ±0.012 for the extraction 
of Vus from Ki 3 decays. 
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1 Introduction 


Chiral perturbation theory [1-3] (CHPT) is the effective theory describing the low-energy ex¬ 
pansion of QCD Green functions. It is a fundamental tool in low-energy hadron phenomenology. 
State-of-the-art calculations involve expansions to NNLO (p®) in external momenta and quark 
masses [4]. Given the large number of low-energy constants (LEGs) appearing in the GHPT 
effective Lagrangian to order p® [5-7], in order to retain predictive power it is highly desirable 
to develop a non-perturbative framework to match the effective low-energy description to QGD 
and to estimate the LEGs. In this work we use the 1/Nc expansion framework together with 
the Minimal Hadronic Ansatz [8-10] to study the matching between low and high energies for 
the (SPP) three-point function of one scalar and two pseudoscalar densities in the chiral limit. 

The purpose of our investigation is twofold. First, we wish to extend the 1/A'"c'-motivated 
matching scheme that has proved successful in other cases [8-16] to the (SPP) three-point 
function. The scheme entails the construction of a hadronic interpolation between the known 
low- and high-momentum regimes, dictated respectively by chiral symmetry and by the QGD 
asymptotic behaviour. The approximations involve the choice of the hadronic content and the 
corresponding set of short-distance constraints to be satished. We truncate the spectrum to 
the lowest-lying resonance multiplet per channel, based on the observation that the low-lying 
hadronic spectrum has the largest impact on the LEGs. This choice dehnes our hadronic ansatz 
as the most general meromorphic function with poles corresponding to Goldstone bosons and 
S, V resonances. Goncerning the short-distance behaviour, the leading power in the OPE for 
the {SPP) Green function displays anomalous scaling, making the explicit matching with a 
meromorphic function problematic. In our analysis we therefore impose that the behaviour of 
our ansatz at large momenta is not worse than the one required by QGD. These constraints 
are then fully compatible with the asymptotic vanishing of those hadronic form factors that 
we need to consider for hxing the relevant LEGs. 

Beyond the general aspects mentioned above, the (SPP) Green function is of special phe¬ 
nomenological interest, as its low-energy behaviour allows to determine the LEGs governing 
SU(3) breaking in Kg^, decays to 0(p®). In turn, decays offer the possibility of a precise 
determination of the GKM mixing parameter Vus, and thus an accurate test of GKM unitarity 
when combined with knowledge of Vud [17]. Our aim is to explore the quantitative implications 
of our matching framework for Kg^ decays and to assess the attendant uncertainty. 

The material in this paper is organized as follows. In Sec. |2l we describe in detail the 
matching procedure for the {SPP) correlator. We give its low- and high-momentum limits, 
the interpolating form and we discuss the implications for the LEGs. In Sec. OJwe then review 
the status of GHPT calculations of Kg^ form factors and the impact of our Endings on the local 
contribution of 0{p^). Finally, in Sec.|3]we summarize our main results and conclusions. 
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2 The (SPP) Green function 

Our starting point is the Fourier transform of the octet (SPP) Green function in massless 
QCD. On account of SU(3) and C invariance it is given by a single scalar function^, 

jdxdyPP^+^^y+^^^{Q\TS^{x)P\y)P%z)\Q) = d^^^Tlspp{p ^, ( 1 ) 

with p + g + r = 0. Bose symmetry implies that the function is symmetric in its second and 
third arguments 

^spp{s,u,t) = Ilspp{s,t,u) . (2) 


2.1 Chiral symmetry 

The low-energy expansion of the function If^pp may be worked out in CHPT: 

^sppis-it,u) = -^-T—^|-Fg -|- 4 L 5 s-|-4(4L8 — L 5 ) (t-|-u) (3) 

tu t 

—8G12 -\- 8(2Gi 2 + G34 -|- Gas) sit + u) — 8(Gi2 -l- G34 — Gas) iP + '^^)j' 

~(4-Bo)^(2Gi2 — 4G31 — 2G34 -|- 2G38 — G94) -|- 0 {p ^), 

where the loop contributions have been discarded. The leading-order contribution to If^pp 
depends exclusively on the pion decay constant Fq and on the quark condensate, Bq = 
— (0|-uu|0)o/Gq in the chiral limit rUu = xrid = = 0. The coefficients of the higher-order 

contributions are written in terms of the coupling constants Lj and Gj [3,6]^. 

2.2 Asymptotic behaviour 

The Operator Product Expansion (OPE) implies the following short-distance behaviour of the 
(SPP) Green function when s,t,u ^ 00 : 


— (t — 

Usppp t, u) ^ 2(0|hu|0)o -^^ , (4) 

stu 

to leading order in inverse powers of momenta, but to zeroth order in the strong coupling 
constant. In fact, the anomalous dimensions of the (pseudo)scalar currents and of the quark 

^The quark currents are normalized according to 

P“(a;) = q{x)ij5X°‘q{x) , S'“(a;) = q{x)X°-q{x) , 

with (A“A^) = furthermore, = i(A“{AbAG)- 

^In the literature, there exist different conventions as to the normalization of the Ci [6,7]. We prefer to 
work with the Ci of mass dimension —2 since they have the canonical large-behaviour. 
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condensate imply nontrivial Wilson coefficients in iD that ensure the correct scale dependence 
and modify the asymptotic behaviour by logarithms. We come back to this important point 
in subsection El 

In the case where only two of the three coordinates x, y, z approach each other, the analysis 
reduces to the OPE of pairs of quark currents and yields in momentum space^ 

^5Pp(p^g^(p + g)^) =0(g“^) , p&xed,q^oo, (5) 

^spp(p^g^(p + g)^) =-8(0|hM|0)o-|^+ 0(p"^) , g£xed,p^cx). (6) 

pZqz 

There is one additional constraint that follows from chiral symmetry: chiral Ward identities 
and pion pole dominance lead to the relation 

lim n5pp(s, t, s) = Uss{s) - npp(s), (7) 

ZJdq 

where we have introduced the scalar two-point functions 1155 and IIpp according to 

tJdxe^P^^-yH0\TS'^{x)S\y)\0) = 6‘^^Uss{p^) , ( 8 ) 

i [dxe^P^^-y\0\TP^{x)P\y)\0) = d'^’^Uppip^) . 


It is well known that the OPE for T3°“ (x) S’^(y) coincides with the one for TP°‘{x)P^{y) in the 
chiral limit up to terms of order {x — y)~^. We thus deduce 


oV nspp(s, t, s) = 0{s 2) , (9) 

^^0 ZJdq 

when s —>• cx), not conflicting with (jHl). 

A different class of asymptotic constraints comes from considering the behaviour of various 
hadronic form factors at high momentum transfer. Form factors may be obtained from II^pp 
by extracting the residues of the appropriate double poles. As an explicit example, the scalar 
form factor of the pion Fg'^{s) is given by the residue of the double pole at t = u = 0, 


fs) = lim 


tu 


t,u^o (2 BqFq)^ 


n 


SPP 


s, t, u) 


( 10 ) 


The asymptotic condition [18] in this case reads 


lim F^^(s) = 0 . (11) 

S—>-oo 

Similar constraints exist for other (transition) form factors. 

^The term given explicitly does in fact involve the two-point function (0|rA“(a;)P*’(?/)|0). However, for 
vanishing quark masses, that correlation function coincides with its short-distance limit, 

(0|TA“{a;)P'’(y)|0) = -2z(5“'’(0|M'u|0)od^Ao(x - y) + 0{mq) . 
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2.3 Model approximation 


In the large-limit the only singularities in n 5 pp are single-meson poles. In this spirit, and 
truncating the spectrum to one resonance per channel, we construct a meromorphic approxi¬ 
mation to the function II^pp with explicit scalar (S) and pseudoscalar (P) poles (besides the 
pion). 


n 


sv 

SPP 


s, t, 




[Ml 


Po + Pi+ P 2 + P 3 + Pi 
- s\[-t][-u][Ml - t][Ml 



( 12 ) 


where the are polynomials of degree n in s, t, u\ 


Pr 


— / , / , ^n—k,k—l,l S t U 


k=0 1=0 


(13) 


The ansatz in Eq. m represents the most general expression for the given particle content 
(hxing the denominator) that does not violate the short-distance behaviour n 5 pp = 0{p~^) 
required by the OPE 0. The normalization is chosen such that Pq = Cqoo = 1 yields the 
correct low-energy limit. Since the function II^pp is symmetric under the interchange of t and 
u we have in addition 

C-kml ^klm ■ ( 1 ^) 

We are thus left with 21 parameters Cioo, • • •, C 022 to be determined. 

It is straightforward to work out the implications of our model for the chiral coupling 
constants that enter in Eq. o. Aside from cqoo = I, the expressions for the LECs involve the 
coefficients in Pi and, in case of the Ci, also those from P 2 : 


Li'p = 


r^sv 

'^12 


r^sp 

'-'34 


r<sp 

'-"38 


FS 


Ml 


+ Cloo 


j^SP ^ fA 
® 16 


1 1 

+ ClOO + c„.„ 


(16) 


1 Tp2 

tSP _ ^0 

'2Ml ^ 8 


C 2 OO 5 


1 

■ 1 1 ■ 


1 

1 

772 

rSV ^0 

2 

[Ml ^ M 2 J 


[Ml 

0 ?" 

® 16 [ 


m|m2 


3o 


200 


Clio + C020 


1 1 
+ 


M| M2 


-SP 


-Cq 

16 


MlMl 


C200 ~ Clio ~ C020 


The coupling constants C 31 and C 94 are not hxed individually. For the combination C 31 -|- IC '94 
one hnds 

1 


1 F2 

risp , ^ r^sp _ r'SP , r'SP ^0 
'-'31 “T ^'-'94 — '-'34 w '-'38 ^2 


K 


+ 2 Co 20 ~ Coil 


(16) 


Before moving on to the asymptotic constraints, let us observe that the relation between 
n^pp and 1155 — IIpp given in Eq. m implies that the combination cioo + coio relevant for Lg 
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is given in terms of quantities that specify the two-point functions. In fact, dehning Cm and 
dm in terms of the one-particle matrix elements of the scalar and pseudoscalar currents [19], 


|(0|^“|5')| = 6'^HV2BoCm , K0|P“|iP')| = 6^^4V2Bodm , 


one derives from Eq. o 


(17) 


11 8 
M| + Mj + ^ 

which implies the well-known result [19] = 1/2 (c^/M| — 

The (SPP) Green function has also been studied in Ref. [14] within a ladder resummation 
inspired hadronic model. Phenomenological consequences such as the LECs of 0(p®) were not 
considered in that approach. 


Ml Ml 


(18) 


2.4 Implementing asymptotic constraints 

The aim of the present investigation is to obtain information on the chiral coupling constants by 
exposing the ansatz in Eq. m to suitable constraints implied by QCD asymptotic behaviour. 
In particular, it has been proven successful [8-16] to match a meromorphic representation 
to the QCD short-distance behaviour, given in this case by Eqs. ©, ®, and ®. Here we 
encounter a problem related to the missing Wilson coefficients in those relations. For instance, 
matching to Eq. would imply that the nine coefficients in P 4 are hxed completely. 


2 C211 — C022 — — 2 c 031 — 


1 

2BIMIMI ’ 


(19) 


with all other coefficients in P 4 being zero. 

The problem alluded to is manifest in Eq. (USD, which would imply that non-vanishing 
coefficients Ckim depend on the running scale of QCD in the same manner as the (quark 
condensate)"^. Including the appropriate Wilson coefficients in Eqs. © and © would repair 
the scale dependence but it would make the Ckim momentum dependent at the same time inval¬ 
idating our ansatz. Moreover, that momentum dependence would involve the strong coupling 
in the non-perturbative regime. It is also obvious that the logarithmic dependence induced by 
nontrivial Wilson coefficients can never be matched by a meromorphic approximation with a 
hnite number of resonances (see, e.g.. Ref. [20]). 

Therefore, we have to hnd an alternative set of criteria to determine the parameters of our 
model ansatz and the relevant LECs. Of course, this discussion is not limited to {SPP) but 
applies to all Green functions with anomalous scaling to leading power in inverse momenta. 

1. While matching onto the short-distance result leads to the problems discussed above, 
we certainly wish to ensure that the short-distance behaviour of our ansatz is not worse than 
what is predicted in QCD. The behaviour n 5 pp = 0(p"^) (up to logarithms) was built into our 
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ansatz from the start. There are, however, nontrivial restrictions that follow from the limits 
in which the momenta are treated asymmetrically: 


IlsPp{p^qMp + qf)=0{q-^), 

p fixed , q —>■ 00 , 

0) 

^5pp(p^ {p + qf) = 0{p -^), 

q fixed , p —>• cx) . 

0) 

Our model (jl2|l reproduces this behaviour, provided the following relations hold: 


C400 + C310 + C220 + Ci30 + Co40 

= 0, 

(20) 

8C4 oo + 3C3 io — Ci30 

= 0, 


6C400 + 3C3 io + C220 

= 0, 


C3IO + C211 + C121 + Co3i 

= 0, 


2 Co 40 + 2 Co 31 -|- C022 

= 0, 


C300 + C210 + C120 + Co30 

= 0 . 



2. More conditions arise from the short-distance behaviour of 1155 ( 5 ) — npp(s) in First, 
the parameters Cm and dm must satisfy [11,21] 


dl 



( 21 ) 


a relation analogous to the hrst Weinberg sum rule [22], In addition, Eq. Q requires that the 
following combinations of coefficients vanish: 


C200 + Clio + C020 = 0 , 
C300 + C210 + C120 + C 030 = 0 , 
C400 + C310 + C220 + Ci 30 + C 040 = 0 . 


( 22 ) 


3. The asymptotic vanishing of the scalar form factor of the pion requires 

Cioo = C200 = C300 = C400 = 0 . ( 23 ) 

When combining the restrictions from the pion scalar form factor in Eq. (OH|l . the scalar/ 
pseudoscalar Weinberg sum rule in Eq. (EH), the OPE conditions in Eq. EOD, and Eq. m , we 
are left with a total of 9 undetermined parameters: 

P2 = ciio[s(t + u) - u^] + contu , ( 24 ) 

P 3 = C2io[s‘^(t + u) - u^] + cnistu + Co2i[t + u]tu + Ci2o[-s(t^ + V?) u^] , 

Pi = C3lo[s^(t + u) — 3s^(f^ -I- M^) -I- 3s{t^ + u^) — (f^ -I- 3tu + u^){t — uY] 

+ C21l[s^ - (^ - uY]tU + Ci2l[s(t + u) - {t- uY]tu . 
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4. Once the conditions from one-pion transition form factors {7i\S\V) (csio = C210 = 0, 
Clio = —MpCi 2 o) and (7r|P|5) (ci2o = 0) are included, the ansatz reduces further to 

P2 = Coiitu , ( 25 ) 

-P3 = [cillS + Co2l{t + u)]tu , 

Pa = [c2ii{s^ - {t - uf) + Ci2i{,s{t + u) - {t - uf)]tu . 

The constraints discussed so far are sufficient to hx the LECs of order and p® within our 
scheme. 

We could £x additional parameters in the polynomials (123) by also considering the various 
form factors of the S and P resonance states. As pointed out in Ref. [14], this procedure in 
general leads to a representation in conflict with the OPE constraints. We do not dwell on this 
issue further since it is of no concern for our purpose of determining the LECs occurring up to 

0(p®). 


2.5 Low-energy constants 

The constraints enumerated above determine the LECs L5, Lg) (^12, (^34 and Cgg in terms of 
resonance masses and couplings. 


• The asymptotic vanishing of the pion scalar form factor hxes L5 and C 12 through cioo = 
C200 = 0- The scalar form factor dehned in Eq. (113 takes the simple form 


= (26) 

s 

The coupling constants describing the dependence on the variable s in the low-energy 
expansion of n5pp, viz. L5 and C 12 in Eq. (ED, are determined by the momentum depen¬ 
dence of the scalar form factor alone. Within the single-scalar resonance approximation 
that implies (see also [23]) 


U = 


F‘2 

4M| 


^12 = 


P2 

np 


8M| 


(27) 


In the vector sector the analogous consideration leads to the well-known predictions 
[12,13,16,24] 


La — 


np 

2M2 


Cg 


C'gp — 


F2 

^0 




(28) 


• Enforcing the correct short-distance behaviour of 1155(5) — npp(s) by virtue of the hrst 
of Eqs. determines the value of Cgg. 


• The asymptotic vanishing of one-pion transition form factors (7r|S'|'P) and (7r|P|5) implies 
Clio = C020 = 0 and thus determines C34. 
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The combination C^i + 1/4(794 remains undetermined, as the coefficient Cqh is not hxed by 
the constraints we have considered. 

Within our framework the relevant LECs of (9(p®) are determined in terms of the scalar 
and pseudoscalar octet masses (M5 and Mp), the pion decay constant Fq and the couplings 
Cm and dm dehned in Eq. (ED- Using also the Weinberg-like sum rule (EH, we obtain 


-sv 


r'SV 

'-'12 


C. 


•sv 

34 


r'SV 

'-^38 


-^0 tSV _ \ ( \ 

4M| ’ 8 “ 2 \Ml Ml,) ’ 

-^0 

8M| ’ 

_LV 

16Af| 2 VjWJ M}j ’ 

Fl ,<{1 _Lf 

16M| 2 VM| j ' 


(29) 


To estimate these couplings, we need numerical values for the input parameters. For the pion 
decay constant we use the physical value Fq = = 92.4 MeV. The coupling dm can be hxed by 

studying the pion scalar form factor away from the chiral limit [9], resulting in dm = Fq/(2^/2) 
(or Cm = -fo/2), which we take as central value. As for the mass parameters, spectroscopy and 
chiral symmetry [25,26] suggest a central value Mp = 1.3 GeV, while Ms is more controversial. 
The analysis of Ref. [26] would suggest Ms = 1.48 GeV for the lightest scalar nonet that 
survives in the large-Vc limit. This result is supported by recent lattice calculations (see for 
example Ref. [27] and references therein). It implies a value of ~ 10“^, in good agreement 
with most recent hts to the (9(p^) LEGs [28]. 



Ci ■ 10^ GeV2 

• (4vr)^F2 

6Ci ■ 10^ GeV^ 

r^sv 

'-'12 

-4.4 

-0.09 

1.6 

r<sv 

'-^34 

6.6 

0.14 

4.7 

r<sv 

'-^38 

2.5 

0.05 

2.2 


Table 1: Numerical values for the LEGs of 0{p^) in GeV“^ and in natural units 
(47r)“^i/7^ for Ms = 1.25 GeV. The last column contains the variations of the 
Ci{p) for Mrj < /i < 1 GeV. More precisely, we display the quantities 6Ci = 
max{|Gi(Mp) — Gj(Mp)|, \Ci{Mp) — Ci{l GeV)|}, using the T[(/i) from ht 10 in Ref. [28]. 

In the estimates below, we use Mp = 1.3 GeV and we vary Ms between 1 and 1.5 GeV. 
The numerical values for Ms = 1.25 GeV are collected in Tabled The above results represent 
the leading term in the large-Vc expansion of the couplings (within our simplihed scheme 














of truncating the spectrum to the lowest-lying resonances). One way to estimate the size 
of sub leading corrections in 1/Nc is to look at the renormalization scale dependence of the 
couplings. This effect is formally higher order in I/Nq and a leading-order estimate is unable 
to provide the scale at which the expressions (|2ni) apply. At O(p^), resonance saturation works 
well for /i = Mp. On this basis, a crude estimate of the uncertainty is given by the variation 
of Oj(/r) for /i between and 1 GeV. 

In this way, we obtain the uncertainties shown in Table ^ Two comments are in order 
here. First of all, it is not the uncertainty of any given LEG that matters^ but the overall 
scale dependence of a measurable quantity, as will be discussed below for the vector form 
factor at t = 0. Secondly, the values in Table ^ are of course very sensitive to the scalar 
resonance mass, while the scale dependence stays the same. This strong dependence on Ms 
seems rather discomforting at hrst sight but, as above, physical observables may exhibit a 
smoother dependence as we will demonstrate in the following section. 

3 SU(3) breaking in decays and V^s 

We now investigate the consequences of our results for the estimate of SU(3) breaking in Ki^ 
form factors. As is well known, decays offer one of the most accurate determinations of 
the GKM element Vus- After the recent re-evaluation of radiative corrections [29,30] and the 
new experimental results [31-34] (see Ref. [35] for a review of the present experimental and 
theoretical status), the main uncertainty in extracting Vus comes from theoretical calculations 
of the vector form factor (0) at zero momentum transfer dehned by 

{7r-{p^)\s'ypU\K^{pK)) = /+°^”(t) {pK+Pn),^ + /5°^'(t) (PK-Pn),, , (30) 

where t = {px —P-kY- Here we are interested in the SU{3) breaking corrections to (0). 

We break up the form factor in terms of its expansion in quark masses: 

f+ ^ (0) = 1 + fpi + fp6 + ... . (31) 

Deviations from unity (the octet symmetry limit) are of second order in SU{3) breaking [36]. 
The hrst correction arises to O(p^) in GHPT: a hnite one-loop contribution [37,38] determines 
fp4 = —0.0227 in terms of TV, Mk and M^r, with essentially no uncertainty. The p® term 
receives contributions from pure two-loop diagrams, one-loop diagrams with insertion of one 
vertex from the effective Lagrangian, and pure tree-level diagrams with two insertions from 
the p^ Lagrangian or one insertion from the p® Lagrangian [39,40]: 

fp‘ = /y'“'”(A>) + + fpSrw . (32) 

Individual components depend on the chiral renormalization scale p, their sum being scale 
independent. Formally speaking, the three contributions scale as 0(1/A^"^), 0{1/Nc), and 

After all, the choice of LECs is basis dependent. 
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0(1), respectively. Although our main concern here is with the other terms have to 

be accounted for in a consistent phenomenological analysis, as infrared logs tend to upset the 
1/Nc counting. Using as reference scale /r = Mp = 0.77 GeV and the Li from £t 10 in Ref. [28], 
one has [40]: 

jyoops^Mp) = 0.0113, (33) 

fL^xioop^M^) = -0.0020 ± 0.0005 . (34) 

Note that we have subtracted the tree-level piece proportional to L 5 x L 5 from the corresponding 
quantity A(0) in Ref. [40]. 

The explicit form for the tree-level contribution is then [40,41] 


f;r{Mp) = 8 


{Ml - Ml 
F3 


{Li{Mp)y 

F3 


- Cl,{Mp) - c;,{Mp) 


(35) 


Upon substituting Ll{Mp) 


Lf’ and Cl^^ziiMp) Gf 2^34 one gets 


fT{Mp) = 


{Ml - 
2 Ml 


Ml 

Ml 


(36) 


where we used dm = F.„/{2\/2). In Fig. [T] we plot as a function of Ms our full result for 
fl^^^{Mp) (solid line) as well as its two components: the L 5 x L 5 term (dashed line) and the 
^12 + ^34 piece (dotted line). The two contributions tend to largely cancel each other, reducing 
the full result to ~ 10 % of each individual term. As a consequence, the ambiguity related to 
Ms is strongly reduced, given the size of the resulting effect. This is simply a consequence 
of treating all tree-level contributions to fp6 on the same footing, as suggested by the l/Nc 
counting. 

In addition to varying Ms in the range 1 GeV < Ms <1.5 GeV, we need to estimate 
the intrinsic uncertainty due to our choice of p = Mp as matching point for /p 6 ®®(Mp). We 
do this by performing the matching for any p G [M^, 1 GeV] and then running the result 
back to /i = Mp via the renormalization group [7]. This way we find 5/p6^®(l/Vc) = ±0.008. 
Accounting also for the uncertainty in Ms, we get 


f^iMp) = - 0.002 ± 0.008 i/Nc ± 0 . 002 Ms ■ (37) 

We now discuss the two main features of our result: 


i. The smallness of the local contribution /* 6 ®®(Mp). 

ii. The size of formally sub leading terms in the 1/Nc expansion: the scale dependence of 
the LEGs and the loop contribution to fp6. 
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1 1.1 1.2 1.3 1.4 1.5 
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-(C12 + C34) 


Ms (GeV) 

Figure 1: We display according to Eq. (p!)|l as a function of Ms for Mp = 1.3 GeV 

(solid line). We also plot the two components according to Eq. the dashed line represents 
the term proportional to L 5 x L 5 , while the dotted line represents the term proportional to 

~(Gi2 + G34). 


The naive expectation for the size of the local contribution to fp6 is /Mg ~ 

5-10“^ GeV^/Mg. An order of magnitude is lost through the factor 1/2(1 —< 0.1 (for 
resonance masses in the range considered) in Eq. . This extra suppression is a consequence 
of imposing, within the particle content of our ansatz, the correct asymptotic behaviour for 
the two- and one-pion form factors (7r|S'|7r), (7r|5|'P) and (7r|P|5). 

We do not have a complete answer to the question whether this suppression persists when 
using a more sophisticated ansatz. However, we have examined the stability of our result in 
two different directions. Omitting the pseudoscalar resonances altogether in our ansatz m 
gives rise to a solution that is equivalent to setting = 0 in Eqs. (Ei. There is a com¬ 
plete destructive interference in this case: the scalar contributions cancel in Eq. implying 
/p6®®(Mp) = 0 instead of Eq. (IHF)|l . On the other hand, adding an additional pseudoscalar nonet 
leads to a straightforward generalization of Eq. (jsni). Instead of 


d 


2 

m 



with 



(38) 


one gets by including a second pseudoscalar multiplet with mass Mpi and coupling 


d 


2 

m 



+ d: 


r/2 



with (ij. + = P ■ 


(39) 


Assuming d'Y < Mg = 1.25 GeV and taking Mp/ in the range 1.5 —>• 2 GeV, \f^l^^{Mp) \ may 
increase by up to 0.002. Thus, additional pseudoscalar multiplets do not modify the general 
result of a small tree-level contribution f^l'^^{Mp). 
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The second issue concerns higher-order corrections in the I/Nq expansion. We have taken 
the variation of the LECs with the renormalization scale as a measure of those corrections. 
The relatively big variation for /i in the range < /^ < 1 GeV is due to large infrared logs 
that typically occur in the scalar sector. This also explains why loop contributions cannot be 
neglected in this case. 

Adding all uncertainties linearly, but omitting the small error of the loop contribution (ED, 
we get hnally® 


fTW 

= -0.002 ± 0.008 ± 0.002Ms p' . 

(40) 

fp^ 

= 0.007 ±0.012, 

(41) 

/f'"(0) 

= 0.984 ±0.012. 

(42) 


Our hnal result (I12D differs from other determinations [38,42,43]. This difference is due to 
the small central value for which appears to be a generic consequence of a few-resonance 
approximation. When combined with the size and sign of the loop contribution of 0{p^) [40], 
our central value for fps in m is positive in contrast to most other estimates, e.g., = 

—0.016 ± 0.008 of Leutwyler and Roos [38]. 

For the purpose of illustration, we use the recent branching ratio measurements of 
KTeV [31], together with their results on form factors where the curvature in (t) 

has been included in the analysis [32]. Using also the recent precise determination of the Kl 
lifetime [33], one hnds [44] (0) ■ |14s| = 0.2166 ± 0.0010. Our result for (0) in 

then implies 

|K,| = 0.2201 ± 0.0027/_,(o) ± O.OOlOe^p . (43) 

The central value is smaller than the one from CKM unitarity, using the most recent value 
of Vud [45]. It is also smaller than what one would obtain using the same experimental input 
and f^°'^ (0) from Refs. [38,42,43]. On the other hand, our result is in better agreement with 
alternative extractions of 14* from r decays [46] and [47]. 

The LECs obtained in this article determine also the deviation from the original Callan- 
Treiman relation [48] 

Act = (Mj - Ml) - ^ , (44) 

involving the scalar form factor fo{t) = f+{t) + tf-{t)/{M]^ — M^). The tree-level contribution 
of 0{p^) is given by [42] 


A tree,p° 

LXct 




Ml) {2 a 


+ Cl,) 


Ml{Ml - Ml 

Ml, 


Ml 

Ml 


(45) 


® Estimates of the uncertainty due to higher-order corrections beyond 0(p®) [40] essentially do not modify 


the final result for /4 ^ (0) when adding the errors in quadrature. 
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With the values of Ms and Mp considered before, is negative and small (a few times 

10“^ in magnitude). Finally, these LECs also provide a prediction for the slope® of the scalar 
form factor. In the representation for Aq given in Ref. [40], the LECs appear again in the 
combination 2 Ci 2 + C 34 . Including the loop contributions [40], we obtain 

Ao = 0.013 ± 0.002 i/Mc ± 0.001 Ms = (13 ± 3) ■ lO’^ , (46) 

in agreement with Aq = (17 ± 4) ■ 10“® from Ref. [37], Aq = (15.7 ± 1.0) ■ 10“^ from Ref. [42] 
and with the value measured by KTeV [32] in decays, Aq = (13.72 ± 1.31) ■ 10“^. 

4 Conclusions 

We have constructed a meromorphic approximation to the {SPP) Green function with pole 
singularities corresponding to Goldstone modes and lowest-lying scalar (S) and pseudoscalar 
(P) resonances. The highlights of our analysis are: 

• We have shown how the polynomial terms in our ansatz can be hxed by imposing the 
correct large-momentum behaviour implied by QCD both off-shell (OPE constraints) 
and on the relevant hadron mass shells (form factor constraints). Because of nontrivial 
Wilson coefficients, one cannot match the OPE constraints exactly with a hnite number 
of resonance poles. We have instead required a large-momentum behaviour for our ansatz 
that is not worse than predicted by QCD. For the particle content used, no inconsistencies 
arise with this set of constraints. 

• This matching procedure has allowed us to determine three of the LECs of 0(p®) (Cu, 
C 34 , Css) appearing in the low-energy expansion of (SPP) in terms of resonance masses 
Ms p and F.^. In particular, C 12 is uniquely determined by requiring the correct behaviour 
of the pion scalar form factor (ttIS'Itt), while C 34 is hxed by the correct scaling of the 
one-pion form factors {n\S\V) and ( 7 r|P|i 5 ). 

• We have estimated the uncertainty of the large-iVc matching procedure by varying the 
chiral renormalization scale at which the matching is performed. While being a sub¬ 
leading effect in the l/Nc counting, the scale ambiguity M,, < /i < 1 GeV gives rise 
to sizable uncertainties. For the same reason, loop contributions must be included in 
phenomenological applications. 

• We have explored the impact of our results on the estimate of the local p® contribution 
to SU{3) breaking in Kis decays. We hnd that the resulting effect is much smaller than 

®The slope Aq is defined by 

/o(t) =/+(0) 1 -I- Ao—h ... . 

-*7r+ 
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the ratio of mass scales (M|- — /Mg would suggest, due to interfering contributions. 
When combined with known loop corrections [40], our estimate leads to fpS = 0.007 ± 
0.012, the mean value being opposite in sign to most other existing calculations. Using 
this input and the most recent experimental results for the branching ratio [31] and 
lifetime for Kl [33], we hnd |14s| = 0.2201±0.0027/_^(o)±0.0010exp- The mean value of our 
result is smaller than the value inferred from CKM unitarity = 0.2265 ± 0.0022, 

using the most recent determination of Vud [45]. If the recent precision measurement of 
the neutron lifetime [49] were confirmed our preferred value for Vus would however be in 
perfect agreement with CKM unitarity. 

• We have considered variations of the hadronic ansatz to investigate the stability of our 

results. The smallness of the tree-level part compared to the loop contribution of 0(p®) 
for (0) appears as generic feature of a few-resonance approximation for the set of 

large-momentum constraints considered. In view of the significant implications for the 
determination of Vus, fhe validity of our approach will be further investigated also for 
other Green functions. 

• Finally, we have also used our results to estimate the deviation from the Callan-Treiman 
relation and to calculate the slope of the scalar form factor. 
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